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Abstract 



We propose a model of a relativistic string formed by a scalar complex 
field, acting as electromagnetic field source. An axiosymmetric solutions 
of the stationary equations for the scalar and electromagnetic fields are 
' found numerically. The mass m is calculated as a function of the charge e 

and the magnetic moment p of the system. The resulting toroidal struc- 
ture is interpreted as an electron because the calculated ratio e 3 /(2mc 2 /i) 
coincides with the fine-structure constant a = e 2 /(he) ~ e 3 /(2m e c 2 /i e ). 
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1 Introduction 



Removing divergences in quantum electrodynamics by renormalization proce- 
dure allows calculating multiloop contributions to electrodynamic quantities 
[1]. As a result, the theory leads to an excellent agreement with experimental 
data. For example, taking into account radiation corrections into the electron's 
magnetic moment leads to an agreement with experiment within 10 significant 
digits. However, the fine-structure constant, which enters all calculations, is 
not calculated within the theory. It is possible that one has to go beyond elec- 
trodynamics in order to calculate it. The possible value 1/(4 tt) of the unificld 
coupling constant in the grand unification schemes at the energies from 10 to 
10 19 GeV is still under discussion (see for example [2]). In addition, calculations 
connecting this constant with the fine-structure constant in the energy interval 
of 17 orders of magnitude seem to be ambiguous [3]. Alternative numerological 
approaches (see for example [4]) give good agreement with experiment but fail 
to make any predictions. 
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In quantum electrodynamics, the point-like electron is surrounded by a cloud 
of electron-positron pairs and is actually a finite-size structure. This fact sug- 
gests that it might be possible to use the finite size structure as a ground state 
approximation in the perturbation theory and thus avoid any divergences in the 
theory. The "combined" models, where two or more material fields are con- 
nected by the vector fields, are well known. For example, an electromagnetic 
field connects electrons with nuclei in atoms, gluon fields connect quarks in nu- 
cleons and mesons [5]. In contrast with this, in a model presented here we have 
only one material field as a source of electromagnetic field. The aim of the paper 
is to investigate the soliton-like states of the system. 

In the previous paper [6], we proposed a classical model of the finite-size 
axiosymmetric particle with charges and currents situated on a disk. These 
charges and currents were assumed to be proportional to the field potentials. 
Kaluza-Klcin nonlincarity of electromagnetic field [7] , [8] ensured the existence 
of stable solutions and finite values of all observables in non-trivial stationary 
states. However, the fine-structure constant calculated within the model turned 
out to be of order one. 

In this paper we propose a new model of extended electron in the classical 
form of a flat ring current. Quantum fluctuations smear this string out to a 
tore, while the charge e and the magnetic moment /i are invariant with respect 
to fluctuations. The basis of the relativistic Lagrangian formalism is a scalar 
complex material field, which acts as a source of the Maxwell electromagnetic 
field. The interacting fields form a stable axiosymmetric soliton-like state. The 
aim of the paper is validation of the proposed relativistic Lagrangian formalism 
for the description of the electron as a finite-size toroidal structure. The system 
of equations derived from this formalism is solved numerically and the energy 
E = mc 2 of the stable state is evaluated. The energy is a function of the model 
parameters e and /i, but the dimcnsionlcss number e 3 / (2 m c 2 /x) is independent 
of them and must coincide with the fine-structure constant a = e 2 /(frc) w 
e 3 / (2m e c 2 fi e ) (if one neglects the electron's anomalous magnetic moment). The 
demonstrated agreement with the experimental value serves as a validation of 
the proposed model. 

2 Lagrangian and parameterization of solutions 

We assume that the electromagnetic field is classical and the Lagrangian has 
the form 




(1) 



F, 



A — A 
dx» v dx v 



The wave function ip is normalized as 





2 



where A is the fundamental length of the model determined later. The interac- 
tion Lagrangian containing only linear terms in potentials is 



Lint — 




ip-z — i>*)A„ d 3 x, 



(3) 



the Lagrangian of the " free" field ip has the form 




x. 



(4) 



We note the parameters of the model are the charge e and the fundamental 
length A. Note also that there are no bilinear in electromagnetic potentials 
terms in Lagrangian components (1, 3, 4), which are required by the gauge 
invariancc. Thus, conservation of electromagnetic current does not follow from 
Lagrange field equations. 

We seek a solution for ip, which is axiosymmetric (around z-axis) and has a 
plane of symmetry z — in the cylindrical coordinate system centered in the 
center of the charge. Conservation of electromagnetic current can be ensured 
by choosing the wave function to be an eigenfunction of the operator ie\ 
and operator ieX/r-^r, such that the phase depends only on t and (j>, and the 
amplitude depends only on r and z: ip = 1/(2tt) 1 ^ 2 iP (r, z) e *0-**A. 

The magnetic moment fj, = /j, e can be obtained by setting A = fj, e /e = 
h/(2m e c) — 1.93 • 10~ 13 m, i.e. 4ir times smaller than the Compton wave 
length of electron. The parameter A can be interpreted as a ring radius of 
the string current at which the current four-vector is light-like. Because the 
fine-structure constant e 3 /(2 to c 2 /i) we are calculating is dimensionless, we set 
e=l,A = l,c=l and calculate the mass to = E. 

Because ip* (r, z) = ip (r, z), the currents components J r and J z are equal to 
zero, and the electromagnetic field has only two non-zero components Ao (r, z) 
and A,), (r,z). Taking all these assumptions into account, calculating the inte- 
grals with respect to <j> and varying the Lagrangian over Ao (r, z), A$ (r, z), and 
ip (r, z), we get the Maxwell equations and the wave function equation: 



V 2 A) (r,z) = -2ip 2 (r,z), 

V 2 yl (r,z) = -2-p 2 (r,z)/r, 
(-l/2X7 2 + U{r,z))iP{r,z) = 0, 



(5) 



(6) 
(7) 




The amplitude of the wave function is normalized as 




(8) 



where dS — rdr dz. 
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Because of the symmetry, the boundary conditions at z = require that all 
functions are even functions of z. The electric field potential is an even function 
of r, and the wave function and the magnetic potential are odd functions of r. 
At large r and z, the wave function should vanish, and the asymptotic behavior 
of Aq and is respectively determined by e and \i. 

The wave function equation (7) has the form of the stationary Schrodinger 
equation with zero self-energy. The magnetic field creates the potential well 
surrounded by the barrier created by the electric field at large distances and 
by the "centrifugal" forces at small distances. The depth of the well, which is 
necessary for existence of a stable solution, can be achieved at a characteristic 
radius of the current that is much smaller than the fundamental length. This 
leads to a small value of the fine-structure constant. 

The wave function decreases at large distances. Therefore, the asymptotic 
behavior of the electromagnetic field potentials are the same as those of the 
point-like particle. This in turn allows calculating the wave function asymptotic 
behavior: 

V> - cxp(-(8 1 / 2 (r^ + z 2 ) 1 / 4 )) r^^/^+z 2 ) 1 ' 2 (l + l.5468/(r 2 +z 2 ) 1/i + -)- (9) 

The relatively slow decay of the wave function (in comparison with Gaussian 
decay exp(— x 2 /2) for the harmonic oscillator wave function) is due to the fast 
decay at the boundaries of the potential well and its smallness in the vicinity 
of the barrier. This fact validates the quasi-classical approximation for the 
ground state. We give a quasi-classical estimate of the toroidal structure radius 
r . The potential well A^/r = 1/(2 7q) and kinetic energy (7r/(2r )) 2 must be 
approximately equal. This implies ro = 2/tt 2 . The energy of the bounded state 
is 

E = J dL/d(d ^*) d a r Sx + J dL/d{d^) d ^ (fx ~L M -L - L mt (10) 

and taking the normalization condition into account, we obtain that 

E = 1 — L M — L — L int . (11) 

For xjj and A satisfying the equations of motion, we have Lo + Li nt = and 
Lint = — 2 * Lm- This results in 

E = l + \ J^ 2 A^/rd 3 x. (12) 

Integrating over (f> leads to the final result for the energy of the bounded state, 

E=l + \ J*P 2 (r,z) A^(r,z)/rdS. (13) 

So, in the quasi-classical approximation we find E = 7r 6 /16 and 1/a = ir 6 /8 w 
120. The fact that such a crude estimate gives a very reasonable result is 
encouraging. 



4 



3 Numerical solution algorithm 



We propose the following change of variables in order to take the non-polynomial 
(exponential) form of the wave function into account effectively: 

r = L r f(p), (14) 

z = L z g(Q, (15) 

where L r and L z are the scale factors, f(p) and g(Q are some logarithmic 
functions which transform infinite quadrant r = O..oo, z = O..00 into the unit 
square p — 0..1, ( = 0..1: 

f(p)= S p+(-ln{l-p t )Y t - i y t /p t -\ (16) 

g{Q = {-ln{l-C)) {v - 1),v /&- 2) . (17) 

After this transformation (at fixed s, t, and v) the wave function is interpolated 
by polynomials, which are easy to integrate and differentiate. The problem then 
split into two parts: 

1) calculating L z at fixed L r /L z , f(p) and g((), which leads to a convergent 
solution and 

2) choosing an "optimal" solution by optimizing L r /L z , f(p) and g((). 
3.1 Calculation of the scale factor 

In order to calculate L z , we note that the nonhomogeneous Maxwell equations 
have nontrivial solution for any wave function in the right hand side, whereas the 
homogeneous Schrodinger equation has only trivial solution with the approxi- 
mate potentials obtained from the approximate solution of Maxwell equations. 
We construct an iterative process by adding nonhomogeneous term proportional 
to ipi to cq. (7) for tpi+i: 

(-1/2 V 2 + Ui (r, z)) V> i+ i (r, z) + U t (r, z) fa (r, z) = 0. (18) 

We start with some initial value of the normalized wave function ip = ipo, 
solve nonhomogeneous Maxwell equations, and then find the next iteration ipi+i 
from eq. (18). The resonance increase of the norm of V'i+i (by 7-8 orders of 
magnitude) indicates that the nonhomogeneous term is small in comparison 
with the analogous term containing Vi+i; an d therefore, the iterations can be 
stopped (after each iteration, the solution is normalized according to eq. (8)). 
In order to calculate the eigenfunction of the wave function equation, we vary 
the scale factor L z (at fixed L r /L z ) until we find the resonance in the norm of 
tp. The energy value is then calculated. 

We put the described procedure to a test by calculating of the odd wave 
function of the first exited state of the 1-d harmonic oscillator. A few first odd 
moments of the probability density coincided with the analytical values within 
a few tenths of percent (due to the interpolation over 8 points). 
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3.2 "Optimal" lattice 

The search of the resonance value of L z was repeated after changing the ratio 
L r jL z or the functions f(p) and g(C). The energy was found to depend (within 
a few percents) on the particular choice of L r /L z , f(p) and g((), suggesting a 
need for an increase of the number of parameters of the wave function. Such an 
approach leads to a prohibitive calculation time. An alternative approach used 
in this paper is to choose the ratio L r /L z and functions f(p) and g(Q in such a 
way that the result achieves its asymptotic value with a relatively small number 
of parameters (of order of thousand) determining the wave function and the 
potentials. For this purpose, we introduce additional requirements on functions 
(14,15), which reduce arbitrariness and increase the accuracy of the result. 

1) We require a balance of contributions to the norm of the wave function 
from large and small distances along r axis and a balance of contributions from 
large distances between r and z axis. This requires that the probability density 
as a function of 2 p — 1,( should be close to axiosymmetric function (i.e., look 
like the up-side down cup) . This is analogous to the requirement that a few first 
moments should be zero (we limited ourselves with the dipolc, quadrupole and 
octopole moments): 

r (2p-l)iP 2 (r,z) dS = 0, (19) 



/< 



/ 



((2p-lf-z 2 )^ 2 (r,z)dS = 0, (20) 

J {2p-lfip 2 {r,z) dS = 0. (21) 

2) We characterize the fields by their values on the lattice with coordinates 
defined by the zeros of the Legendre polynomials in order to use the Gaussian 
quadratures numerical integration method. We can also estimate the accuracy 
of obtained solution by calculating the probability of the total square deviation 
of the full solution from the solution ^4 interpolated over each fourth point of 
the lattice: 

5Pi = j (V (r, z) - V> 4 (r, z)f ^ 2 (r, z) dS. (22) 

The minimum of this deviation, i.e. the good approximation of the solution by 
the low order polynomial, suggests a high precision of numerical integration and 
differentiation. 

3) We can impose one more condition on the ratio L r /L z : we assume that 
the characteristic volume L 2 L z is less sensitive to numerical errors and require 
that its derivative with respect to the ratio L r /L z (for fixed f(p), g(() and 
resonance Lz) is zero: 

° ■L 2 r (L r /L z )L z (L r /L z )=0. (23) 



d(L r /L z ) 

We should stress that conditions (19. .23) are not physical equations but serve 
as a criterion for the most accurate solution at fixed lattice size. We impose 



G 



these conditions one by one and calculate the ratio L r / ' L z and the energy H 
for fixed f(p) and g(C)- The stability of results (minimum of relative dispersion 
of L r /L z ) with respect to the choice of conditions 19.. 23 is considered as an 
evidence of a good choice for functions f(p) and g{(). 

The most time-consuming part of the calculation is solving the system of 
linear equations for the field values on the lattice sites. The calculation time 
grows as the cube of the number of lattice sites. It also grows with the increasing 
precision of numerical differentiation. Therefore we differentiate the wave func- 
tion using its interpolation over 9 or 11 points and electromagnetic potentials 
over 5 points because they are smoother than the wave function. Numerical in- 
tegration using Gaussian quadratures is exact for polynomials of order of twice 
number of points on the lattice. In this work we use lattices 20 x 10, 24 x 12 
and 28 x 14 along the respective p and Q axis. 

4 Numerical results 

We use transformations (14, 15) and seek a resonance value of the scale pa- 
rameter L z , which leads to a convergent solution of Maxwell equations (5, 6) 
and the wave equation (7) for fixed functions (16,17) and different values of the 
ratio L r /L z . Using five additional criteria (19. .23) for the optimal lattice, we 
calculate five different values of the ratio L r /L z . A small scatter of these values 
6L/L = S(L r /L z )/(L r /L z ) is considered as a criterion for the optimal choice of 
functions (16, 17). The table below lists the results for three lattices of different 
size where the minimum value of 6L/L is below 3%. 

First two columns specify the number of lattice sites along p and £ respec- 
tively and the number of points used for calculating derivatives of function tp. 
Third column contains mean square deviation of smoothed solution from the 
exact one(22). Next column lists relative deviation of the ratio L r /L z from its 
mean value for 5 different optimal lattice criteria. Fifth column contains main 
results of the paper: the product of twice the energy and magnetic moment 
divided by the charge cubed. 



Lattice sites 


Num. diff. 


SPi 


SL/L 


2 to c 2 /i/e 3 


r 


V z 2 


20 x 10 


9 pt 


0.1576 


2.1 % 


134.13± 0.14 


0.1362 


0.0744 


24 x 12 


11 pt 


0.0144 


1.4% 


136. 55± 0.20 


0.1363 


0.0735 


28 x 14 


11 pt 


0.0073 


0.3 % 


136. 82± 0.05 


0.1363 


0.0718 



The experimental value of the quantity 2m e c 2 p e /e 3 is 137.19 (it is slightly 
higher than 1/a = 137.03 because of the anomalous magnetic moment of the 
electron). So, the calculations on the 28 x 14 lattice agree with the experimental 
value within few tenths of percent, and the discrepancy decreases with increasing 
the lattice size. As expected, the characteristic size of the toroidal structure is 
almost an order of magnitude smaller than the fundamental length (see two last 
columns of the table, which list the ratios of the average radius and height of 
the toroidal structure to the fundamental length) and two orders of magnitude 
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smaller than the Compton length of electron. Its typical values are of order 
2 • 10~ 14 m. The correction to the second bounded state of the toroidal structure 
in non-uniform proton field is 7 orders of magnitude smaller than the second 
level energy of a hydrogen atom and 10 times smaller than the Lamb shift. Thus, 
the contribution from the finite size of the totoidal structure is comparable with 
the radiation corrections, at least at small energies. 

5 Conclusions 

We propose a relativistic Lagrangian formalism, which lead to a system of equa- 
tions for the complex scalar and electromagnetic fields. Assuming axial and 
mirror symmetries and specific angular and time dependence of the wave func- 
tion, we obtain a system of three equations for the stationary amplitude of the 
wave function and two electromagnetic potentials as functions of two spatial 
variables. We construct an iteration procedure for solving these equations on a 
lattice and the algorithm for selecting an optimal lattice. The numerical solution 
of this system allow calculating the energy (or the mass) of the soliton-like state. 
We interpret the obtained toroidal structure as an electron because the value 
of dimensionless constant e 3 /(2mc 2 /i) coincides with the experimental value 
e 3 /(2 m e c 2 /j, e ) w e 2 / (he) = a within few tenths of percent. Further increase of 
numerical precision is desirable. 

One of possible confirmations of the proposed Lagrangian formalism and 
the analogous axial symmetry might be the construction of another soliton-like 
solution with larger mass and smaller magnetic moment that can be interpreted 
as ^i-meson. 

The correspondence principle requires constructing the effective Lagrangian 
of toroidal structure, investigating it's properties, including gauge invariance, 
and solving the well known electrodynamic problems on its basis. One could 
try to find electron stationary states in external fields, in particular, to calculate 
bound states of toroidal electron in the field of the proton (hydrogen atom) . 

I am grateful to I.V. Tyutin, B.L. Voronov and A.E. Likhtman for useful 
discussions. 
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